Sign reversal of the order parameter in s-wave superconductors 
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^ ! Abstract 
l> ■ 

We show that in a superconductor where two or more bands cross the Fermi 

o ■ 

l/^ ' level it is possible, in the framework of the conventional (s-wave) BCS theory, 

o\ ' 

' that the sign of the superconducting gap is different on the different sheets of 

^ ' the Fermi surface. At least one of two conditions has to be satisfied: (1) The 

■ 

' interband pairing interaction is weaker than the Coulomb pseudopotential, 

o : 

^ ' while the intraband one is stronger, or (2) there is strong interband scattering 

> '■ 

' by magnetic impurities. In the case of YBa2Cu307 we shall argue that the 

^' 

' first condition is possibly satisfied, and the second one very likely satisfied. In 

many aspects such a sign-reversal s— wave superconductor is similar to a d- 
wave superconductor, and thus demands revising recent experiments aimed to 
distinguish between the s— and d— wave superconductivity in this compound. 
1. Recently, a number of experiments probing the relative phase of the order parameter 
A on different parts of the Fermi surface in the superconducting YBa2Cu307 have been 
reported [T|-§[ • Some of them [0,^0 seem to indicate the conventional pairing state with A 
having the same sign over the whole Fermi surface, while others suggest that A changes in 
sign, as consistent, e.g., with the ci-pairing. The question of the symmetry of the supercon- 
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ducting state, and thus of the interpretation of these experiments, is of crucial importance for 
distinguishing between the conventional mechanism for superconductivity and more exotic 
mechanisms, or among the unconventional theories themselves. 

It is generally believed that the s-pairing is inconsistent with sign reversal of A. This is 
not true. A simple counterexample is the case of two concentric Fermi spheres, which have 
gaps (order parameters) of opposite signs. Such a state has pure s— symmetry. The two- 
dimensional analog are two coaxial Fermi cylinders. A Josephson contact between two such 
superconductors, or with a conventional superconductor, may show an unusual behavior, 
sometimes similar to the c?— pairing case. 

In this paper we shall show under which conditions the situation similar to the examples 
above, which we shall call the Interband Sign Reversal of the Gap (ISRG), can be realized, 
and we will also argue that these conditions are not at all exotic but are likely to be realized 
in YBa2Cu307. We shall also discuss briefly how ISRG can manifest itself in Josephson 
tunneling, and we shall make a link to the existing experiments. 

2. The extension of the BCS theory for two or more superconducting bands was first 



worked out by Suhl, Matthias, and Walker M and independently by Moskalenko [|TU|, and 



later elaborated on by many. It was realized |TT| that the fact that several bands cross the 



Fermi level is not sufficient to have considerable many-band effects in superconductivity. 
Only when the bands in question have a very different physical origin, can a substantial 
effect appear. 

This is the case in many high-T^ cuprates. In particular, YBa2Cu307 is known to have 
four sheets of the Fermi surface, all four having a different physical origin |T^: One is 



formed by the chain pda (seen by positron annihilation), another is an apical oxygen band 
(seen in de Haas- van Alphen experiments), and the last two are bonding and antibonding 
combinations of the two pda plane bands (seen by angular- resolved photoemission) . Basing 
on the richness of the band structure of YBa2Cu307, several groups pointed out that at 
least the two-band |T^, or probably the whole four-band [ p!^P^ , picture should be used to 



describe superconductivity in this system. Various experiments have been interpreted as 



indicating two or more different superconducting gaps. 



We shall now remind the basic equations of the multiband BCS theory P,p!T|: The 



Hamiltonian has the following form: 




where e^^fc is the kinetic energy in the z-th band, c*;,^ and Ci^ka are corresponding creation 



and annihilation operator, and gij is the averaged pairing potential. 

The order parameter A on the z-th sheet of the Fermi surface is given by the equation 



if the cut-off frequency cu/j is assumed to be the same for all sheets. is defined in the 
usual way by the effective coupling constant, \og{2'y*Lj£)/nTc) = l/Xejf, 7* — 1.78. The 
effective coupling constant Ae// in this case is simply the maximal eigenvalue Amax of the 
matrix Aij = QijNj, where Nj is the density of states at the Fermi level (per spin) in the j-th 
band. Aij plays the role of the coupling constant A in the one-band BCS theory. Note that 
conventional (isotropic) A is also defined in terms of Aij: A = -^ijNi/N = Yl,i\Ni/N ^ 
where the mass renormalization for the i-th band is Aj = ^ij) and N = J2i ^i- Obviously 
Ae// > A, which means thatdue to larger variational freedom Tc in the multiband theory is 
always larger than in the one-band theory. The two are equal in isotropic case, i.e. when Qij 
does not depend on An instructive example of the opposite case is the two-band model 
with All = A22 = A > 0, A12 = A21 = -A. Then A = 0, while X^jf = 2A. Note that the 
last value is the same as when A12 = A21 = A. The physical reason is that although there is 
no solution of Eq.|l| with Ai > 0, A2 > 0, there is an obvious solution with Ai = — A2 7^ 0. 
Near the solution of Eq.|l| is A2/A1 = (Ae// — Aii)/Ai2, demonstrating directly that the 
sign reversal of the order parameter, A2/A1, takes place when nondiagonal matrix elements 
A12 and A21 are negative. The fact that conventional BCS theory (Eq.|^) allows for the 
ISRG solution, has never, to our knowledge, been mentioned in the extensive literature 
existing on multiband superconductivity [Q.One can easily check that Eq.|l] may have a 




(1) 
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superconducting solution even for all Qij < 0, i.e., when no attractive interaction is present in 
the system. The condition for that is \gi2\ > {\gi\Nf + \g2\N2) /2NiN2- This is similar to the 
well-known fact that in a system with repulsion the superconductivity with higher angular 
momenta (p, d) is possible, because of the sign reversal of the order parameter. The main 
difference is that in the example above the symmetry of the superconducting state is the same 
as of the normal state. Below we shall demonstrate that even a fully attractive interaction 
Qij > can lead to the sign reversal if (a) interband pairing interaction is weaker than 
Coulomb preudopotential, (b) there is strong interband scattering by magnetic impurities. 

3. If (7's are electron-phonon pairing potentials, then Eq.|l] should be corrected for a 
Coulomb repulsion, which can be readily done [jll| by substituting gij — > gij — U*j ~ 



gij — U*, where the effective Coulomb repulsion U* is logarithmically renormalized in the 
same way as in one-band superconductivity theory {U* is assumed to be independent on 
A direct consequence of that is that if the interband electron-phonon coupling is 
weak, the situation with a negative gap, g^j — U* < 0, can easily be realized because of 
the interband repulsion. We illustrate that by numerical calculations presented in Fig.|^. In 
this calculations the following parameters had been used: gi2 = g22 = 0, A^i = 4:N2, and 
gii = Ni'^ so that to have A = 1. This choice of parameters corresponds to the plane and 
chain bands in YBa2Cu307, as discussed below in Section 5. Several facts draw attention: 
First, in this model Tc decreases with the increase of fi* = U*N substantially slower than 
in a one-band case when — * when /i* — > A. Second, the order parameter induced in 
the second band ("chains") is always negative; Its absolute value reaches maximum when 
|Ai| = IA2I, I.e., at U* = gppNp/2{Np - N,). 

4. Following the standard way of including the impurity scattering in the BCS theory 
||17|| , one writes the equations for the renormalized frequency iD„ and order parameter A^ 
{n is the Matsubara index), which completely define the superconductive properties of the 
system (see, e.g. Ref. |T8[| ): 



A.n = A. + E^(7.--7^-) 
Ai = TTTj2AijA,JQj,n. (2) 



Here u;„ = (2?2 + 1)tcT, Qi^n = y^ln + ^i,ny lij is the scattering rate from band i into band 
j due to nonmagnetic impurities, and Yij is the same for magnetic impurities. Near Eqs.H 
can be solved analytically [|19|. For two bands, in the linear in 7, 7s approximation, the 
solution redices again to Eqs. |1|, with the effective coupling matrix A: 



^ 27n + 1I2 + 712 1I2 - 712 ^ 



■A. 



y 721 - 721 27I2 + 7 + 721 J 

When all A's are equal, the standard Abrikosov-Gorkov result is recovered: 5\ ~ — 7rA^(7j^ + 
7i2 + 721 + 722) /8^co- The main point of the AG theory [1^ is that 7"* enters equations for uj 



and A with opposite signs. That is why the magnetic impurities appear to be pair-breakers, 
and the non-magnetic ones not. The above solution shows that in the multiband case of 
Eqs.0 this argument works only for the intraband non- magnetic scattering (7jj drop out), 
while all other scattering rates are, in principle, pair-breaking. 

An interesting special case is Ai2,A2i ^ Aii,A22- Then in the effective A matrix non- 
diagonal elements Alj^^i 7^ j) = A^- + nAiiAjj{jij — 7|^)/8Tco can become negative, if 
•y-j is sufficiently large. As discussed above, this situation will lead to ISRG. In order to 
demonstrate this effect quantitatively, we solved the EqsJ^ in the Eliashberg approximation 
numerically, using the following parameters: An = 1,A22 = 0.5, A12 = 0.025, A21 = 0.1. 
This choice is again inspired by the situation in YBa2Cu307:The ratio of the densities of 
states in the bonding and antibonding bands in YBa2Cu307, according to the band struc- 
ture calculations, is about 2.5 but it is likely the calculations underestimate this value 
(see discussion below). Correspondingly, we used 721 = 4712,721 = 47i2. The results for the 
low-temperature regime, T << Tc, are shown in Fig.^ In accord with the condition derived 
above, when the difference 712 — 712 becomes larger than some critical value (in this case, 
0.042'7rTc), the second gap changes sign. In other words, when attractive interband coupling 



is relatively weak and the magnetic interband scattering is strong the system will choose to 
have two gaps of the opposite signs, losing in pairing energy, but avoiding the pair-breaking 
due to interband scattering. 

5. Let us consider now two cases relevant for YBa2Cu307_5, and their applications for 
the Josephson effect. 

(1) It is believed by many (e.g., Ref. [|l^) that the chain electrons would not be su- 
perconducting, or only weakly superconducting, if not for the "proximity effect" from the 
planes. In our language, this means that Ac — /ic = (fi'cc — U*)Nc ~ 0, where c stands for 
chains. Then the sign of the gap, induced in the chains, will be determined by the sign of 
Ope — U*; in a quite likely case of Qcc < U*, ISRG between the chain and the plane bands 
takes place. 

(2) One can also look for the ISRG between the bonding (6) and antibonding (a) com- 
binations of the two plane bands. According to the calculations ||2^ and experiment ||21| , it 



is the a band which has Van Hove singularities near the Fermi level. In the calculations the 
singularities are bifurcated, which makes the density of states in the antibonding band 2.5 
times larger than in the bonding band, and is exremely sensitive to the warping of Cu02 
planes, thus resulting in strong electron-phonon interaction. Experimentally, the singulari- 
ties are even closer to the Fermi level than in the calculations, and are extendended towards 
F-point. If, as it is often claimed, this singularity plays a crucial role in superconductivity, 
then the a band is the superconducting one, and the superconductivity in the b band is 
induced. Consequently, one has the situation similar to the above- described "j>-c" scheme. 

(3) Furthermore, the ISRG due to magnetic impurities may also be relevant for 
YBa2Cu307. Let us assume that the main magnetic scatterers are antiferromagnetic (AF) 
spin fluctuations on the plane Cu sites. For the moment we assume these fluctuations to be 
static (see, however, the discussion below). Inelastic neutron scattering studies |2^ show that 
the AF correlations between the planes survive even in the fully oxygenated samples, where 
the intraplanar correlations are virtually non-existent (correlation length ^/a = 0.84 ±0.04). 
This is in direct contradiction with the popular assessment that the intraplane AF correla- 
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tions are more important than those between the planes. To understand the consequences of 
this fact we shall again consider the bonding and the antibonding band. The former is even 
with respect to the z — > —z reflection, and the latter is odd. The standard Hamiltonian for 
the magnetic scattering is 

^l,- = - EE < - R)SRfT|j/5 >, 

R a/3 

where Sr, is the spin of the impurity at point R, and are spin indices. In case of 
two antiferromagntically correlated impurities in the two planes, and (j75| must be of 
different parity to render non-zero HI-. This means that only r^^ is non-zero. 

In the previous paragraph we considered static impurities. In this context "static" means 
that the characteristic frequency of the AF fluctuations Tiioaf ~ ttTc =25 meV. It is not clear 
yet how large loaf is in YBa2Cu307_2,.. Detailed calculations will have to include the proper 
frequency dependence of the magnetic susceptibility in the same way as it is done with the 
phonons in the Eliashberg theory. It is obvious, however, that the qualitative conclusions 
will not change. 

6. Let us discuss the consequences for the Josephson effect separately for each of the 
above cases. Rather than trying to explain the contradictory experimental results reported 
so far we shall indicate some qualitative predictions of our model. 

The supercurrent density through a grain boundary in a two-band superconductor can be 
written as Jg = Jl^ sin ^^■'^where J^^is the Josephson critical current density correspond- 
ing to tunneling between the bands iand j,and ^^■'is the gauge-invariant phase difference of 



order parameters Aj and A,-. In the simplest case of the Ginzburg-Landau regime [OB 



where Rij = {h/e'^){pF/27T?i)'^/{D) is the tunneling resistance per unit area, (D) is the 
angle-averaged transparency of the barrier and pp = min{pFi,pFj). Following Geshkenbein 
and Larkin P^, we obtain immediately that if the order parameters have different signs, 



sign{Ai) = —sign{Aj), then in the stationary case (J^ = 0) the finite phase difference 



appears, (pij = n. This is similar to the "vr— contact" considered by Bulaevskii et ah p5 
but in our case it is due to the sign reversal of the order parameter in different bands. 

Generally, the total critical current, J*°*, depends on orientation of the boundary relative 
to crystallographic axes because of the angular dependence of Rij. It can become negative for 
certain directions when the contribution due to interband tunneling prevails. The condition 
is + J2i > J]i + J22 for the HTS/HTS junction and J^^ > J}^ for the HTS/LTS one. To 
some extent this effect is similar to that considered by Sigrist and Rice for the c/— pairing, 
but some of our predictions differ qualitatively, as discussed below. 

a) Let us first consider the c — p scenario. For a HTS/LTS junction the tunnel resistance 
Ruid) depends strongly on the angle 6 relative to the b axis, namely Ruid) has a sharp 
minimum at 6 = due to the strong angle dependence of a barrier transparency (D) for a 
tunneling process. Moreover, according to band structure calculations kinetic energy of 



carriers along the chains is larger than that in plains, thus leading to larger (D) values. As 
a result, J*°*(6')<0 for small 6, whereas for all other angles J*°*(6')>0. Therefore an intrinsic 
vr- phase shift will occur in this case between tunneling along a-and along 6-directions. Then 
a dc SQUID with junctions on a and b faces of a crystal, will show a $o/2 shift of a field 
dependence Ic{H). This effect was observed in Refs. [0,^0 and attributed to the d^2_y2 
pairing state. Another consequence is a shift of a Fraunhofer pattern for a single junction 
formed on the corner of a crystal, because Jc changes sign along the junction , as was 
discussed in Refs. [|l],^ for d - pairing. Evidently, the same effect will take place in the 
considered case. 

We shall also mention that nonzero Josephson current, observed for c-axis tunneling in 
Pb/insvAaXoi/ Yi^xPi"xBa2Cu30j tunnel junctions in Ref. ^ is contradictory to the d^i-^yi 
symmetry. Indeed, nonzero Josephson current was predicted theoretically for c-axis contact 
between s-wave and d-wave superconductors in Ref. [^], but only in the second order in 
boundary transparency (D) , therefore this model can not explain the rather large values of 
IcRn products of the order of 1 meV observed in Ref. On the other hand, this observation 
is consistent with the suggested c — p scenario. The reason is that, contrary to the case of 
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dx2_y2 symmetry, an average order parameter in the ab plane is nonzero. 

Interesting consequences appear for HTS/HTS (grain-boundary) junctions. As follows 
from the above arguments, if 6* = in only one of the grains, then the grain boundary 
is a TT-contact, otherwise it is a conventional one. Consider a closed contour crossing 
grain boundaries. The flux quantization condition in zero external field reads n$o = LPj + 
X]m=i,Af('^o/27r)0-J*'', where L is a self-inductance of a ring and 0lj*^is a phase difference 
across m-th junction. Then it follows immediately, that if a contour crosses an odd number 
of 6 = junctions, a spontaneous magnetisation of a ring with half- integer flux quantum 
will occur, and when it crosses an even number of 6 = junctions the flux quantum will be 
integer. 

Spontaneous magnetization with half-integer flux quantum in a three-junction ring and 
with integer flux quantum in a two-junction ring was demonstrated recently for YBa2Cu307 
in Ref 0. In this experiment all grain-boundaries were of 6* = type. Thus, the results 
are in agreement both with our proposal and with the cZ-wave scenario discussed by Sigrist 
and Rice. To distinguish between these two explanations the measurements for different 
grain orientations are necessary. At the same time, the absence of angular dependence of Jc 
observed in Ref. for a number of different grain-boundary orientations in YBCO does not 
contradict our scheme. Indeed, in all six grain boundaries have had 6^0 which results 
in Jl°\0) -Jl^{e) ~ const. 

Another interesting phenomenon observed first in Bi-based HTS [^and more recently 



in YBCO is the paramagnetic Meissner effect ("Wohlleben effect"). The explanation 



was proposed in Refs. [^,^,^] in terms of intrinsic vr-junctions between weakly coupled 
superconducting grains, giving rise to spontaneous orbital currents in arbitrary directions. 
An external magnetic field will align those spontaneous current loops and can produce a 
net positive magnetization. Therefore, existence of vr-junctions between at least some of 
superconducting grains is a key point for the Wohlleben effect. As discussed above, such 
intrinsic vr— junctions may exist in the considered two band superconductor with interband 
gap sign reversal, thus leading to the possibility of the Wohlleben effect. 
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b) For our second scenario a sign of the total critical current in an HTS/LTS and an 
H TS/HTS iunction in any given direction depends crucially on a relation between the cur- 
rent components J^^ and J^^, i.e. on corresponding tunnel resistances. From symmetry 
considerations, there exists no fundamental reason for a sign change of J^"^ — J]^ vs 9 in 
the ah plane. As a result, the intrinsic phase shift between the two bands can not be de- 
tected by Josephson experiments, similar to the discussed above c-p scenario. In all junction 
geometries such an ISRG superconductor will behave like a conventional s-wave one. 

However, a possibility of the Wohlleben effect still exists: Each given contour would 
include an even number of 7r-contacts, which in an ideal case would compensate one another. 
But in reality the tunnel resistances of these contacts, depending of the local state of each 
grain boundary, will be different, so that the compensation would become incomplete. 

We note, that this scenario may be relevant not only to YBCOhni also to all double-plane 
materials, like i?i-based compounds. 

7. The main goals of the current paper were to demonstrate the possibility of the 
existence of an s-wave superconductor with the sign-reversal of the gap, to point out some 
factors which favor such a state in YBa2Cu^0'j, and to emphasize that some experiments 
interpreted as unambiguous evidences for the (i-wave pairing can, in fact, be explained by 
the suggested ISRG s-pairing. 

In addition to that, we would like to outline some implication of our analysis to a model of 
spin-fluctuation induced superconductivity (Ref. |Q and references therein). In this model, 
AF fluctuations are dynamic, and serve as the intermediate bosons to give superconductivity. 
Only one plane is considered, and, much in the same spirit as in our analysis, the order 
parameters A(k)and A(k + Q), where Q is the AF vector (^, |), are of the opposite signs, 
which in the case of the YBa2Cu^0j Fermi surface leads to the — y^symmetry for A(k). 
Apparently, if one considers two AF coupled planes, and two bands, a and 6, then only in the 



a-h channel does a non-zero pairing potential appear. This is similar to the observation ||31 
that in YBCO the gerade (with respect to z ^ —z) phonons contribute to the intraband 
{a-a and h-h) coupling only, and the ungerade phonons to the interband coupling only. 
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Even without solving the corresponding equations, one can immediately predict the results: 
Since in a bilayer there is no problem having gaps of the same sign on a given sheet of the 
Fermi surface, an s-wave solution must exist, with Aq and Af, having opposite signs. Direct 
numerical calculations show indeed that Montoux-Pines model for a bilayer has a stronger 
instability in s-channel with ISRG than in d-channel Another interesting point is that 
the intraband phonon pairing and the interband spin-fluctuation pairing can coexist and 
even help each other. 

Acknowledgements. The authors are thankfull to O.V.Dolgov, G.M.Eliashberg, 
D.I.Khomskii, A.I.Liechtenstein, E.G.Maksimov and M.R.Trunin for many useful and stim- 
ulating discussions. A.A.G. acknowledges partial support by Russian State Program HTSC 
under project N 93-194. 



11 



REFERENCES 

[1] D.A. WoUman et alPhys. Rev. Lett, 71, 2134 (1993). 

[2] P. Chaudhari and S.-Y. Lin, Phys. Rev. Lett, 72, 1084 (1994). 

[3] D.A. Brawner and N.R. Ott, preprint. 

[4] I. Iguchi and Z. Wen, Phys.Rev.B,49, 12 388 (1994). 

[5] A.G.Sun et al, Phys. Rev. Lett. 72, 2267 (1994). 

[6] C.C.Tsuei et aL, preprint. 

[7] J.Buan et al., Phys. Rev. Lett. 72, 2632 (1994). 
[8] S.Riedhng et aL, Phys. Rev. B 49, 13 283 (1994). 

[9] H. Suhl, B.T. Matthias, and L.R. Walker, Phys. Rev. Lett., 3, 552 (1959). 
[10] V.A. Moskalenko, Fiz. Met. i Met. 4, 503 (1959). 

[11] W.H. Butler and P.B. Allen, in: Superconductivity in d- and f-metals, ed. by 
D.H.Douglass (Plenum, N.Y.,1976). 

[12] O.K.Andersen et al. Physica C (1991) 

[13] V.Z.Kresin and S.Wolf, Phys. Rev. B, 41, 4278 (1990). 

[14] I.I. Mazin et al, Physica C209, 125 (1993). 

[15] L.Genzel et al Z. Phys. B, 90, 3 (1993). 

[16] For a mathematically similar case of a layer system with a spatially varying gap, this fact 
has been pointed out by N.Lazarides, T.Schneider and M.P.Sorensen, Physica C 210, 
228 (1993); For the superconductivity in a two-band Hubbard model due to magnetic 
interaction, ISRG was found by N. Bulut, D.J. Scalapino, and R.T. Scalettar, Phys. 
Rev. B 45, 5577 (1992). 

12 



[17] A.A. Abrikosov and L.P. Gor'kov, JETP 12, 1243 (1961). 

[18] P. Entel et al Z. Phys. B, 363-368 (1975). 

[19] I.I.Mazin, A.A.Golubov, unpublished. 

[20] O.K.Andersen et al, Phys. Rev. B49, 4145 (1994). 

[21] A.A. Abrikosov, J.C. Campuzano, K. Gofron, Physica C214, 73 (1993). 

[22] J. Rossat-Mignod et al in Frontiers in Solid State Sciences, ed. L.C. Gupta and M.S. 
Multani (World Scientific, Singapore, 1993), vol. 1, p. 265. 

[23] V.Ambegaokar and A.Baratoff, Phys.Rev.Lett. 10, 486 (1963). 

[24] V.B. Geshkenbein and A.I. Larkin, JETP Lett. 43, 395 (1986). 

[25] L.N.Bulaevskii, V.V.Kuzh and A.A.Sobyanin, JETP Lett. 25, 290 (1977). 

[26] M. Sigrist and T.M. Rice, J. Phys. Soc. Jpn. 61, 4283 (1992). 

[27] Y.Tanaka, Phys. Rev. Lett. 72, 3871 (1994). 

[28] W.Braunisch et al., Phys. Rev. Lett. 68, 1908 (1992). 

[29] D.I.Khomskn, J. Low Temp. 95, 205 (1994). 

[30] P. Monthoux and D. Pines, Phys. Rev. B49, 4261 (1994). 

[31] O.K. Andersen, unpubhshed. 

[32] A.I.Liechteinstein, I.I.Mazin, and O.K.Andersen, to be published. 



13 



FIGURES 

FIG. 1. Critical temperature and superconducting gaps in a system with induced superconduc- 
tivity in the second band, as a function of Coulmb pseudopotential 

FIG. 2. Superconducting gaps at T = Tc/2 in a two-band system with interband scattering on 
magetic and non-magnetic impurities. Note the straight hne corresponding to A2 = 
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